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Abstract 

In this paper we analyze crowd turbulence from both classical and quantum perspective. 
We analyze various crowd waves and collisions using crowd macroscopic wave function. In 
particular, we will show that nonlinear Schrodinger (NLS) equation is fundamental for quan- 
tum turbulence, while its closed-form solutions include shock-waves, solitons and rogue waves, 
as well as planar de Broglie's waves. We start by modeling various crowd flows using classical 
fluid dynamics, based on Navier-Stokes equations. Then, we model turbulent crowd flows 
using quantum turbulence in Bose-Einstein condensation, based on modified NLS equation. 

Keywords: Crowd behavior dynamics, classical and quantum turbulence, shock waves, soli- 
tons and rogue waves 
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1 Introduction 



Massive crowd movements can be today precisely observed from satellites. All that one can see 
is physical movement of the crowd. Therefore, all involved psychology of individual crowd agents 
(and their groups within the crowd): cognitive, motivational and emotional, as well as its global 
sociology, is only a non-transparent input (a hidden initial switch) to the fully observable crowd 
physics [1] [31 131 H]. About a decade ago, D. Helbing discovered a phenomenon called crowd 
turbulence (see El 13 ISl IH] ) , depicting crowd disasters caused by the panic stampede that can 
occur at high pedestrian densities and which is a serious concern during various disasters (bushfires, 
tornados, earthquakes), as well as mass events like soccer championship games or annual pilgrimage 
in Mecca. 

The adaptive, wave-form, nonlinear and stochastic crowd dynamics has been modeled using (an 
adaptive form of) nonlinear Schrddinger equatioT^ (NLS), also called Gross-Pitaevskii equatior^ 
(GP), defining the time-dependent complex-valued macroscopic wave function tp = ■0(a;,t), whose 
absolute square |'0(a;,i)P represents the crowd density function. In natural quantum units (h = 
1, TO = 1), our NLS equation reads: 

idtr^^-\d,^^ + V\4^\'i>, (i = V^; with a.^ = |^), (1) 

where V = V{w,x) denotes the adaptive heat potential (trained by either Hebbian or Levenberg- 
Marquardt learning). Physically, the NLS equation ([T]) describes a nonlinear wave in a quantum 
matter (such as Bose-Einstein condensates). 

For general crowd simulation, we recently proposed two NLS-based approaches in [1], each 
representing a quantum neural network [10 : 

Weak coupling approach: A looser (and more abstract) but higher-dimensional approach con- 
sisting of an n— dimensional set of NLS equations: 

ydt^,^~\d^.i^, + V\^,\^i^,, (i-l,...,n), (2) 

which self-organize in a common adaptive potential V — V{w,x). Here, the squared am- 
plitude IVjp is the condensate density. The potential V{w,x) includes synaptic weights w^, 

^Thc most important case of nonlinear Sehrodinger equation is the cubic NLS 

with the cubic nonlinearity IV'PV'- The sign -|- in the cubic NLS represents defocusing NLS, while the — sign 
represents focusing NLS. This extraordinarily rich nonlinear PDE represents a fully integrable Hamiltonian system 
which is traditionally studied on Euclidean domains R" (but other domains, like circle, torus or hypersphere, are also 
studied) - and allows for a 'zoo' of various wave- like solutions (to be analyzed later). Terms sub-critical, critical, 
and super-critical are frequently used to denote a significant transition in the behavior of a particular equation 
with respect to a specified regularity class (or conserved quantity). Typically, sub-critical equations behave in an 
approximately linear manner, supercritical equation behave in a highly nonlinear manner, and critical equations 
are very finely balanced between the two. Occasionally one also discusses the sub— criticality, criticality, or super- 
criticality of regularities with respect to other symmetries than scaling, such as Galilean invariance or Lorentz 
invariance. For survey of recent advances in nonlinear wave equations based on their criticality, see |12l . 

^NLS or GP is also similar in form to the Ginzburg-Landau equation, a mathematical theory used to model 
superconductivity. 
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Figure 1: Simulating random-mixing behavior of four crowds in a confined environment, 
which iteratively update according to the standard Hebbian rule: 

n 

Wi ^ -Wi + c\^\g^\il!\, V{w,x) = y^Wjffi, 

where c is the learning rate parameter and gi are Gaussian kernel functions with means Xi 
and standard deviations normalized to unity. The system ^ was numerically solved using 
the Method of Lines (combined with the fast adaptive Runge-Kutta-Fehlberg integrator with 
Cash-Karp accelerator) . Its solutior|f| represents time evolution in the complex plane C of n 

•^Crowd simulations were based on the following data: 

1. Target function used for this case is: / = 2sin(207rt). An infinitesimal tracking function used is: h = 
-{PDF ■ dx/x). 

2. Gaussian kernel functions are defined as: gi = exp[—{v — rrii)^], where v = f — h, while rrii are random 
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cooperative groups, each consisting of m agents of S'i?(2)-kinematic type0 In this approach, 
each individual hue (or kinematic trajectory), defines a velocity controller for a single agent. 
The total number of agents, as well as number of groups, is limited only by the available 
computation power. 

Strong coupling approach: A pair of strongly-coupled NLS equations with Hebbian learning: 
BLUE : i9t^B = - yI'^rI'^--^b + ^IV'bIVb, 
RED : = -^.\^^\'^d,^^^^ + v\cj,^\^^, (3) 



HEBB : w, = - + ch|V'bI5#rL ^ = I] 



Here, ob, foR, cr are parameters related to Red, Blue and Hebb, equations respectively. This 
is a bidirectional, spatiotemporal, complex-valued associative memory machine, generalizing 
Lanchester and Lotka-Volterra predator-prey dynamical systems, as well as Hopfield, Kosko 
and Grossberg models of neural networks (see, e.g. [IT]). 

In this paper we will analyze crowd turbulence (from both classical and quantum perspective) 
as well as various crowd waves and collisions. In particular, we will show that NLS equation ^ is 
fundamental for quantum turbulence, while its closed-form solutions include shock-waves, solitons 
and rogue waves. Firstly, we will model various crowd flows using classical fluid dynamics, based 
on Navier-Stokes equations. Then, we will model turbulent crowd flows using quantum turbulence 
in Bose-Einstein condensation, based on modified NLS equation. 



2 Classical Approach to Crowd Turbulence 

In this section we model various crowd fiows using models from classical fiuid dynamics, based on 
Navier-Stokes partial differential equations (PDEs). 



2.1 Classical Turbulence and Crowd Flows 

Turbulence has long been one of the great mysteries in nature, with discussion dating back to the 
era of Leonardo da Vinci. He observed the turbulent fiow of water and drew pictures showing 



Gaussian means. 

3. Potential field update is given by the scalar product: Vi+i = Vi + Wi ■ gi. 

4. The complex plane C was embedded in the 3D graphics environment (both urban and bush) with 3D collision 
dynamics. 

■^The Euclidean motion group SE(2) = SO{2) X R is a set of all 3 X 3— matrices of the form: 

cos 9 sin 9 x 
— sin 9 cos 9 y 
1 



including both rigid translations (i.e., Cartesian x, j/— coordinates) and rotation matrix 
clidoan plane ss C (see [TsllTe] '). 



cos 9 
— sin 6 



sm y 
cos 9 



in Eu- 



4 



that turbulence has a structure comprised of vortices of different sizes (Figure [2]). After Leonardo, 
turbulence has been intensely studied in a number of fields, but it is still far from completely 
understood. This is primarily because turbulence is a strongly nonlinear dynamical phenomenon. 




Figure 2: Sketch of turbulence by Leonardo da Vinci. 



Turbulent flow is a fluid flow regime characterized by low momentum diffusion, high momentum 
convection, and rapid variation of pressure and velocity in space and time. Flow that is not 
turbulent is called laminar flow. Also, recall that the Reynolds number Re characterizes whether 
flow conditions lead to laminar or turbulent flow. The structure of turbulent flow was first described 
by A. Kolmogorov. Consider the flow of water over a simple smooth object, such as a sphere. At 
very low speeds the flow is laminar, i.e., the flow is locally smooth (though it may involve vortices 
on a large scale). As the speed increases, at some point the transition is made to turbulent (or, 
chaotic) flow. In turbulent flow, unsteady iJorticej^l appear on many scales and interact with each 

^Vortex can be any circular or rotary flow that possesses vorticity. Vortex represents a spiral whirling motion 
(i.e., a spinning turbulent flow) with closed streamlines. The shape of media or mass rotating rapidly around a 
center forms a vortex. It is a flow involving rotation about an arbitrary axis and can be described by the vector 
curl operator. 

In the atmospheric sciences, vorticity is a property that characterizes large-scale rotation of air masses. Since the 
atmospheric circulation is nearly horizontal, the 3D vorticity is nearly vertical, and it is common to use the vertical 
component as a scalar vorticity. 

A vortex can be seen in the spiraling motion of air or liquid around a center of rotation. Circular current of 
water of conflicting tides form vortex shapes. Turbulent flow makes many vortices. A good example of a vortex is 
the atmospheric phenomenon of a whirlwind or a tornado. This whirling air mass mostly takes the form of a helix, 
column, or spiral. Tornadoes develop from severe thunderstorms, usually spawned from squall lines and supercell 
thunderstorms, though they sometimes happen as a result of a hurricane. (A hurricane is a much larger, swirling 
body of clouds produced by evaporating warm ocean water and influenced by the Earth's rotation. In particular, 
polar vortex is a persistent, large-scale cyclone centered near the Earth's poles, in the middle and upper troposphere 
and the stratosphere. Similar, but far greater, vortices are also seen on other planets, such as the permanent Great 
Red Spot on Jupiter and the intermittent Great Dark Spot on Neptune.) Another example is a meso-vortex on 
the scale of a few miles (smaller than a hurricane but larger than a tornado). On a much smaller scale, a vortex 
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other. Drag due to boundary layer skin friction increases. The structure and location of boundary 
layer separation often changes, sometimes resulting in a reduction of overall drag. 

Applied to crowd dynamics, vorticity u) — a;(x, i) is defined as the circulation per unit area at 
a point in the crowd flow field, that is as the curl of the the crowd flow velocity: a; = V x u. ft is 
a vector quantity, whose direction is approximately along the axis of the swirl. The movement of 
a crowd flow can be said to be vortical if the fluid moves around in a circle, or in a helix, or if it 
tends to spin around some axis. Such motion can also be called solenoidal. 

Because laminar-turbulent transition in crowd dynamics is governed by Reynolds number, the 
same transition occurs if the size of the crowd is gradually increased, or the viscosity of the crowd 
is decreased, or if the density of the crowd is increased. 

In particular, in a turbulent crowd flow, there is a range of scales of the crowd flow motions, 
called eddies. A single packet of crowd flow moving with a bulk velocity is called an 'eddy'. The 
size of the largest scales (eddies) are set by the overall geometry of the crowd flow^ 

Such turbulent crowd flow shows characteristic statistical behavior (compare with [SHI HI])- For 
simplicity, we will assume a steady state of fully developed turbulence of an incompressible classical 
crowd flow. Energy is injected into the crowd flow at a rate e in an energy-containing range. In 
an inertial range, this crowd energy is transferred to smaller length scales without dissipation. In 
this range, the crowd is locally homogeneous and isotropic, leading to energy spectral statistics 
described by the Kolmogorov law^ 

E{k) = C e'"^ k-'>'\ (4) 

The crowd energy spectrum E{k) is deflned by _E = J E{k) dk, where E is the kinetic energy of the 
crowd per unit mass and k is the crowd wave-number from the Fourier transform of the velocity 
fleld (compare with tha last subsection on quantum crowd waves). The spectrum of equation ([4]) 
is derived by assuming that E(k) is locally determined only by the crowd energy flux e and by 
k. The crowd energy transferred to smaller scales is dissipated at the Kolmogorov wave-number 
kx = {e/v'^Y^^ in an energy-dissipative range via the viscosity of the crowd flow at a dissipation 
rate £ in equation ([¥]), which is equal to the crowd energy flux $ in the inertial range. The 
Kolmogorov constant C is a dimensionless parameter of order unity. The inertial range is thought 
to be sustained by a self-similar Richardson cascade in which large crowd eddies break up into 
smaller eddies through crowd vortex reconnections. 

In order for two crowd flows to be similar they must have the same geometry and equal Reynolds 
numbers. When comparing crowd flow behavior at homologous points in a crowd model and a full- 
is usually formed as water goes down a drain, as in a sink or a toilet. This occurs in water as the revolving mass 
forms a whirlpool. (A whirlpool is a swirling body of water produced by ocean tides or by a hole underneath the 
vortex, where water drains out, as in a bathtub.) This whirlpool is caused by water flowing out of a small opening 
in the bottom of a basin or reservoir. This swirling flow structure within a region of fluid flow opens downward from 
the water surface. In the hydrodynamic interpretation of the behavior of electromagnetic fields, the acceleration of 
electric fluid in a particular direction creates a positive vortex of magnetic fluid. This in turn creates around itself 
a corresponding negative vortex of electric fluid. 

^For comparison, in an industrial smoke-stack, the largest scales of fluid motion are as big as the diameter of the 
stack itself. The size of the smallest scales is set by Re. As Re increases, smaller and smaller scales of the flow are 
visible. In the smoke-stack, the smoke may appear to have many very small bumps or eddies, in addition to large 
bulky eddies. In this sense, Re is an indicator of the range of scales in the flow. The higher the Reynolds number, 
the greater the range of scales. 

^The Kolmogorov spectrum has been conflrmed experimentally and numerically in fluid turbulence at high 
Reynolds numbers. 
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scale crowd flow, we have Re* — Re, where quantities marked with * concern the flow around 
the crowd model and the other the real crowd flow. 

2.2 Navier-Stokes Crowd Fluids 

Fluid dynamicists believe that Navier-Stokes PDEs accurately describe turbulence (see, e.g. |49]). 
Therefore, we can assume that viscous crowd flows evolve according to nonlinear Navier-Stokes 
equationflj 

li + u- Vu + Vp/p = i/Au + f, (5) 

where u = u(x, t) is the 3D velocity of a crowd flow , u = dtu is the 3D acceleration of a crowd 
flow , p = p(x, t) is the crowd pressure field, while f = f (x, t) is the external nonlinear energy source 
to the crowd, while p, v are the crowd flow density and viscosity coefficient, respectively. Such a 
crowd flow can be characterized by the ratio of the second term on the left-hand side of equation 
([5]), u • Vu, referred to as the crowd inertial term, and the second term on the right-hand side, 
v/^u, that we call the crowd viscous term. This ratio deflnes the Reynolds numherQ Re — vD/v, 
where v and D are a characteristic velocity and length scale, respectively. When v increases and 
the Reynolds number Re exceeds a critical value, the crowd changes from a laminar to a turbulent 
state, in which the crowd flow is complicated and contains eddies0 To simplify the problem, we 
can impose to f the so-called Reynolds condition, (f • u) = e, where e is the average rate of energy 
injection. 

In mechanical Lie algebra terms (see Appendix) , the Navier-Stokes PDE ((S]) can be written: 

a; = — [u, a;] + curl f + vAus. 
On the other hand, the Euler equation for ideal crowd fiows, 

u + u- Vu + Vp/p = 0, (6) 

* Navier-Stokes equations, named after C.L. Navier and G.G. Stokes, are a set of PDEs tliat describe the motion 
of liquids and gases, based on the fact that changes in momentum of the particles of a fluid are the product of 
changes in pressure and dissipative viscous forces acting inside the fluid. These viscous forces originate in molecular 
interactions and dictate how viscous a fluid is, so the Navier-Stokes PDEs represent a dynamical statement of the 
balance of forces acting at any given region of the fluid. They describe the physics of a large number of phenomena 
of academic and economic interest (they are useful to model weather, ocean currents, water flow in a pipe, motion 
of stars inside a galaxy, flow around an airfoil (wing); they are also used in the design of aircraft and cars, the study 
of blood flow, the design of power stations, the analysis of the effects of pollution, etc). 

^Reynold's number Re is the most important dimensionless number in fluid dynamics and provides a criterion 
for determining dynamical similarity. Where two similar objects in perhaps different fluids with possibly different 
flow— rates have similar fluid flow around them, they are said to be dynamically similar. Re is the ratio of inertial 
forces to viscous forces and is used for determining whether a flow will be laminar or turbulent. Laminar flow occurs 
at low Reynolds numbers, where viscous forces are dominant, and is characterized by smooth, constant fluid motion, 
while turbulent flow, on the other hand, occurs at high Res and is dominated by inertial forces, producing random 
eddies, vortices and other flow fluctuations. The transition between laminar and turbulent flow is often indicated 
by a critical Reynold's number (Recrit), which depends on the exact flow configuration and must be determined 
experimentally. Within a certain range around this point there is a region of gradual transition where the flow is 
neither fully laminar nor fully turbulent, and predictions of fluid behavior can be difficult. 

^•^The first demonstration of the existence of an unstable recurrent pattern in a turbulent hydrodynamic flow 
was performed in IKKOll . using the full numerical simulation, a 15,422— dimensional discretization of the 3D Plane 
Couette turbulence at the Reynold's number Re = 400. The authors found an important unstable spatiotemporally— 
periodic solution, a single unstable recurrent pattern. 
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reads in Lie algebra terms, 

d> — —[u,uj], u) — curl u. 

Equation ([SJ is related to the Navier-Stokes PDE ([5]) in the same way as the classical Euler 
equation of a rigid body (with a fixed point, see Appendix), 

■k — TV X <jj, 

is associated to a more general equation, involving friction and external angular momentum j l3U14j 

Tr = 7rxu; + F — vk, (7) 

the 'friction operator' v is symmetric and positive definite. The distributed mass force f , which 
appeared in the Navier-Stokes equation ([SJ, is similar to the external angular momentum F, and 
it is the origin of the crowd flow motion. The viscous friction v/S.\i is analogous to the term ~vk 
in ([7]) slowing the rigid body motion. The similarity becomes especially noticeable if one rewrites 
the equations in components in the eigenbasis of the friction operator. 

For example, for the Navier-Stokes equation with periodic boundary conditions one can expand 
the crowd vorticity field and the force f into the ordinary Fourier series. The equations in both of 
the cases have the form of Galerkin approximation: 

Xi — ^ ^ aijkXjX}^ -f- ^ ^ fi ViXi- (8) 

j,k i 

The first term corresponds to the Euler equation ^ and describes the inertial crowd motion. It 
follows from the properties of ^ that the divergence of this term is equal to zero. Furthermore, 
the Euler equation of an ideal crowd flow in any dimension, as well as that of a rigid body, has a 
quadratic positive definite first integral, the kinetic energy. Therefore, for f = u = 0, the vector 
field on the right-hand side of ((H]) is tangent to certain ellipsoids centered at the crowd origin. This 
implies that during the crowd evolution defined by this equation there is neither growth nor decay 
of solutions. 

The term —ViXi in corresponding to the crowd friction, dominates over the constant 'crowd 
pumping' f when considered sufficiently far away from the crowd origin. Hence, in that remote 
crowd region, the crowd motion is directed towards the origin, and an infinite growth of solutions 
is impossible. Also, since the 'crowd pumping' f pushes a phase point out of any neighborhood 
of the origin, while the friction returns it from a distance, crowd motion in the system of a rigid 
body (O approaches an intermediate regime-attractor. For instance, this crowd attractor can be 
a stable stagnation point or a periodic crowd motion, while for sufficiently high dimension of the 
phase space it can appear to be a chaotic motion sensitive to the initial conditions. 

Here we recall that chaos theory, of which turbulence is the most extreme form, started in 
1963, when E. Lorenz from MIT took the Navier-Stokes PDEs (O and reduced them into three 
first-order coupled nonlinear ODEs, to demonstrate the idea of sensitive dependence upon initial 
conditions and associated chaotic behavior. The 3D phase-portrait of the Lorenz system Q shows 
the celebrated ^Lorenz mask\ a special type of fractal attractor^^ 

^'^The Lorenz reduced system of nonlinear ODEs 

X = a{y — x), y = bx — y — xz, z = xy — cz, (9) 
where x, y and z are dynamical variables, constituting the 3D phase-space of the Lorenz flow; and a, b and c 
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If the crowd friction (or viscosity) coefficient v is higli enough, then the crowd attractor wih 
necessarily be a stable equilibrium position. While the parameter v decreases (i.e., the reciprocal 
parameter, the Reynolds number Re = l/v, increases), bifurcations of the crowd equilibrium are 
possible, and the crowd attractor can become a periodic motion and later a totally 'stochastic' 
oneHl 

2.3 Isovorticial 2D Crowd Flows 

Crowd 2D flow differs sharply from crowd 3D flow. According to V. Arnold, in the realm of fluid 
dynamics, the essence of this difference is contained in the difference in the geometries of the orbits 
of the co-adjoint representation (see Appendix) in the two and 3D cases [131 The character of 
the first, inertia, term in the Galerkin approximation ([S]) (of the Navier-Stokes PDEs (O) changes 
drastically in the passage from 2D crowd flow flows to three- (or higher-) dimensional ones. The 
reason lies in the distinctions among the geometries of the coadjoint orbits of the corresponding 
diffeomorphism groups. In other words, in the 2D case the orbits are in some sense closed and 
behave like a family of level sets of a function. In the 3D case the orbits are more complicated; in 
particular, they are unbounded (and perhaps dense). The orbits of the coadjoint representation 
of the group of diffeomorphisms of a 3D Riemannian manifold can be described in the following 
way. Let vi and be two vector fields of velocities of a non-compressible crowd flow in the region 

are the parameters of the system. Originally, Lorenz used this model to describe the unpredictable behavior of 
the weather, where x is the rate of convective overturning (convection is the process by which heat is transferred 
by a moving fluid), y is the horizontal temperature overturning, and z is the vertical temperature overturning; the 
parameters are: a = P— proportional to the Prandtl number (ratio of the fluid viscosity of a substance to its thermal 
conductivity, usually set at 10), b = ij— proportional to the Rayleigh number (difference in temperature between the 
top and bottom of the system, usually set at 28), and c = K—a number proportional to the physical proportions of 
the region under consideration (width to height ratio of the box which holds the system, usually set at 8/3). The 
Lorenz system ^ has the properties: (i) symmetry: (x, y, z) — >■ {—x, —y, z) for all values of the parameters, and (ii) 
the 2— axis (x = y = 0) is invariant (i.e., all trajectories that start on it also end on it). 

Today, it is well— known that the Lorenz model is a paradigm for low-dimensional chaos in dynamical systems 
and this model or its modifications are widely investigated in connection with modelling purposes in meteorology, 
hydrodynamics, laser physics, superconductivity, electronics, oil industry, chemical and biological kinetics, etc. 

The Lorenz mask (3D chaotic attractor) has the following characteristics: (i) trajectory does not intersect itself 
in three dimensions; (ii) trajectory is not periodic or transient; (iii) general form of the shape does not depend on 
initial conditions; and (iv) exact sequence of loops is very sensitive to the initial conditions. 

^■^The hypothesis that this mechanism is responsible for the phenomenon of turbulenization of a fluid motion 
for large Reynolds numbers has been suggested by many authors. In particular, to normalize the attractor, A.N. 
Kolmogorov suggested in 1965 considering the 'pumping' proportional to the same small parameter u as viscosity, 
and he formulated the following two conjectures |14) : 

1. The weak conjecture: The ma^ximum of the dimensions of minimal attractors (attractor is called a minimal 
attractor if it does not contain smaller attractors) in the phase space of the Navier-Stokes PDEs JSj (as well 
as of their Galerkin approximations l[8)) grows along with the Reynolds number Re = 1/1/. 

2. The strong conjecture: Not only maximum, but also the minimum of the dimensions of the minimal attractors 
mentioned above increases with the Reynolds number Re. As far as we know, both of these hypotheses still 
remain open. 

For the Lorenz system, the role of energy is played by a nonhomogeneous quadratic function. The instability in 
the Lorenz model is apparently stronger than in the Kolmogorov one. One can check how the motion along the 
Lorenz strange attractor sensitively depends on the initial conditions, while for the Kolmogorov model it remains a 
conjecture. It is proven only that a stationary flow indeed loses stability as the Reynolds number Re increases. 

In 1970 Ruelle and Takens formulated the conjecture that turbulence is the appearance of attractors with sensitive 
dependence of motion on the initial conditions along them in the phase space of the Navier-Stokes equation | 17| . 
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D. We say that the fields vi and V2 are isovorticial if there is volume-preserving diffeomorphism 
g : D ^ D which carries every closed contour 7 in D to a new contour such that the circulation 
of the first field along the original contour is equal to the circulation of the second field along the 
new contour: 



Vl — QJ V2 

7 "'97 

The image of an orbit of the co-adjoint representation in a Lie algebra is the set of vector 
fields isovorticial to the given vector field. We have the following law of conservation of circulation 
[T51 [Tl] : The circulation of a field of velocities of an ideal crowd fiow over a closed fiow contour 
does not change when the contour is carried by the flow to a new position. 

Now, for simplicity, we will assume that the region D of the crowd flow is 2D and oriented. The 
metric and orientation give a symplectic structure on D] the field of crowd velocities has divergence 
zero and is therefore Hamiltonian. Therefore, this vector field is given by a Hamiltonian function 
(many- valued, in general, if the region D is not simply-connected). The Hamiltonian function of 
a field of crowd velocities is called the stream- function, and is denoted by ip. Thus we have: 

V = /grad -0, 

where / is the operator of clockwise rotation by 90°. 

The stream function of the commutator of two crowd vector fields turns out to be the Jacobian 
(or the Poisson bracket of Hamiltonian formalism) of the crowd stream functions of the original 
vector fields: 

The vorticity (or curl) of a 2D crowd field of velocities is the scalar function r such that the 
integral around any oriented crowd region ct in of the product of r with the oriented area element 
is equal to the circulation of the field of crowd velocities around the boundary of a: 

rdS — (p V. 

The crowd vorticity can be now computed in terms of the crowd stream function as: 

r = —A-tp. 

In the simply-connected 2D case, isovorticity of crowd vector fields vi and V2 means that the 
functions ri and r2 (the vorticities of these fields) are carried to one another under a suitable 
volume-preserving diffeomorphism. Therefore, if two vector fields are in the image of the same 
orbit of the co-adjoint representation, then a whole series of functionals are equal. For example, 
the integrals of all powers k of the crowd vorticity are: 

r'ldS^ [ r^dS. 

D Jd 

In particular, Euler's equations of motion of a 2D ideal crowd flow: 

dfV + vVv = — Vp; div V = 0, 
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have an infinite collection of first integrals. For example, such a first integral is the integral of any 
power k of the crowd vorticity of the field of crowd velocities: 



where A denotes the antisymmetric wedge product. 

Following V. Arnold [131 [H] , we obtain in this way the following assertions regarding stability 
of planar stationary crowd flows: 

1 . A stationary flow of an ideal crowd flow is distinguished from all crowd flows isovorticial to it 
by the fact that it is a conditional extremum (or critical point) of the crowd kinetic energy. 

2. If (i) the indicated critical point is actually an extremum, i.e., a local conditional maximum 
or minimum, (ii) it satisfies certain (generally satisfied) regularity conditions, and (iii) the 
extremum is non-degenerate (the second differential is positive- or negative-definite), then 
the stationary crowd fiow is stable (i.e., is a Lyapunov stable equilibrium position of Euler's 
equation) . 

3. The formula for the second differential of the crowd kinetic energy, on the tangent space to 
the manifold of crowd fields which are isovorticial to a given one, has the following form in 
the 2D case. Let D he a region in the Euclidean plane with cartesian coordinates x and 
y. Consider a stationary crowd flow with stream function ip = ip{x,y). Then we have the 
quadratic crowd energy form cPH, given by 



where 5v is the variation of the field of crowd velocities (i.e., a vector of the tangent space 
indicated above), and 5r = curl 5v. We remark here that for a crowd stationary flow, the 
gradient vectors of the crowd stream function and its Laplacian are collinear. Therefore 
the ratio AipfVAtp makes sense. Furthermore, in a neighborhood of every point where the 
gradient of the crowd vorticity is not zero, the crowd stream function is a function of the 
vorticity function. 

The above assertions lead to the conclusion that the positive or negative definitenes^ of the 
quadratic crowd energy form cPH is a, sufficient condition for stability of the stationary crowd fiow 
under consideration. The analogous proposition in the linearized fluid dynamics problem is called 
Rayleigh's theorem. 

Now, consider N crowd vortices with velocity circulations fc^, (i = 1, N) around them in the 
Euclidean plane K^. Then the crowd vorticity at any moment will be concentrated at N points, and 
the crowd circulations at each of them will remain constant forever. It is convenient to write the 
evolution of crowd vortices as a dynamical system in the crowd configuration space for the A'^- vortex 

^^A definite bilinear form is a bilinear form B{x,x) over some vector space V such that the associated quadratic 
form Q{x) = B{x, x) is definite, that is, has a real value with the same sign (positive or negative) for all non-zero 
X. According to that sign, B is called positive definite or negative definite. If Q{x) takes both positive and negative 
values, the bilinear form B{x, x) is called indefinite. If B(x, x) > for all x, then B is said to be positive semidefinite. 
Negative semidefinite bilinear forms are defined similarly. 
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system, the space with coordinates Zi = {xi, yi) and symplectic structure kidyif\dxi. Then 
the crowd vortex evolution in M?-^ will be given by the following Kirchhoff-Hamiltonian system 

3 Quantum Approach to Crowd Turbulence 

In this section we model turbulent crowd flows using models of quantum turbulence in Bose- 
Einstein condensation, based on modified nonlinear Schrodinger equation. Here, we want to go 
beyond classical turbulence as described in the previous section using Navier-Stokes equations. 
Essentially, we want to achieve two things here: (i) to provide a cleaner (more controllable and 
repetitive) simulation environment for crowd turbulence; and (ii) the ability to include into this 
environment a variety of nonlinear waves (e.g., shock-waves, solitons, breathers, rogue waves). 
From the physical perspective, the so-called 'Bose-Einstein condensate' will be our macroscopic 
quantum crowd superfluid model. 

3.1 Quantum Turbulence 

Quantum turbulence was discovered in superfluid helium (^He) in the 1950s, but the field moved 
in a new direction starting around the mid 1990s (see [HI [20]). Briefly, quantum turbulence (see 
Figure [3]) is comprised of quantized vortices that are definite topological defects arising from the 
order parameter appearing in Bose-Einstein condensation (BEC)o Hence quantum turbulence 
is expected to yield a simpler model of turbulence than does classical turbulence based on the 
Navier-Stokes PDE (P. 

Bose-Einstein condensation is often considered to be a macroscopic quantum phenomenon (see, 
e.g. [21]). This is because bosons occupy the same single-particle ground state below a critical 
temperature through Bose-Einstein condensation to form a macroscopic wave function (order pa- 
rameter) extending over the entire system. As a direct result of the formation of a macroscopic 
wave function, quantized vortices appear in the Bose-condensed system. A quantized vortesS^ is a 
vortex of inviscid superflow, and any rotational motion of a superfluid is sustained by quantized 
vortices. A quantized vortex is stable and well-defined topological defect, very different from clas- 
sical vortices in a conventional fluid. Hydrodynamics dominated by quantized vortices is called 
quantum superfluid dynamics, and turbulence comprised of quantized vortices is known as quantum 
turbulence (QT) [HJHD]- 

Liquid ^He enters a superfluid state below the A point (T\ — 2.17 K) with Bose-Einstein con- 
densation of the ''He atoms0 The A transition is closely related to the Bose-Einstein condensation 

Bose-Einstein condensate (BEC) is a state of matter of a dilute gas of weakly interacting bosons confined 
in an external potential and cooled to temperatures very near absolute zero. It is the most common example of 
quantum media. 

^^The studies of quantized vortices originally began in 1950s using superfluid *He, and much theoretical, numerical, 
and experimental effort has been devoted to the field. Superfluid ^He, discovered in 1972, presented a system with 
a variety of quantized vortices characteristic of p-wave superfluids. 

^^The characteristic phenomena of superfluidity were experimentally discovered in the 1930s by L. Kapitza 1221 
and Allen et al.\23\. The hydrodynamics of superfluid helium are well described by the two-fluid model proposed by 
Landau |24 | and Tisza 1251 . According to the two-fluid model, the system consists of an inviscid superfluid (density 
p^) and a viscous normal fluid (density p^) with two independent velocity fields Vx and Vn. The mixing ratio of the 
two fiuids depends on the temperature. As the temperature is reduced below the A point, the ratio of the superfluid 
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Figure 3: Quantum turbulence (QT) in Bose-Einstein condensation (BEC). 



of ^He atoms, as first proposed by i23- The Bose-condensed system exhibits the macroscopic wave- 
function V'(x, i) — i)|e'^^'''*^ as an order parameter. The superfluid velocity field is given by 
Vx = {h/m)V9, with boson mass m, representing the potential flow. Since the macroscopic wave 
function should be single-valued for the space coordinate x, the circulation F = ^ v • c?£ for an 
arbitrary closed loop in the fluid is quantized by the quantum k = h/m. A vortex with quantized 
circulation is called a quantized vortex. Any rotational motion of a superfluid is sustained only by 
quantized vortices 

The idea of quantized circulation was first proposed by F. Onsager, for a series of annular rings in 
a rotating superfluid ^28i . R. Feynman considered that a vortex in a superfluid could take the form 
of a vortex filament, with the quantized circulation k and a core of atomic dimension [29jF^ Vinen 

component increases, and the fluid becomes entirely superfluid below about 1 K. The two-fluid model successfully 
explained the phenomena of superfluidity, while it was known in 1940s that superfluidity breaks down when it flows 
fast |26| and this phenomenon was not explained through the two-fluid model. This was later found to be caused 
by turbulence of the superfluid component due to random motion of quantized vortices. 

^'^A quantized vortex is a topological defect characteristic of a Bose— Einstein condensate, and is different from a 
vortex in a classical viscous fluid. First, the circulation is quantized, which is contrary to a classical vortex that can 
have any circulation value. Second, a quantized vortex is a vortex of inviscid superflow. Thus, it cannot decay by 
the viscous diffusion of vorticity that occurs in a classical fluid. Third, the core of a quantized vortex is very thin, 
of the order of the coherence length, which is only a few angstroms in superfluid ^He. Because the vortex core is 
very thin and does not decay by diffusion, it is always possible to identify the position of a quantized vortex in the 
fluid. These properties make a quantized vortex more stable and definite than a classical vortex. 

^^Early experimental studies on superfluid hydrodynamics focused primarily on thermal counterflow. The flow is 
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confirmed Feynman's findings experimentally, by showing that the dissipation arises from mutual 
friction between vortices and the normal flow [32l |33l [SH [35] . Vinen also succeeded in observing 
quantized circulation using vibrating wires in rotating superfluid ^He [3^. Subsequently, many 
experimental studies have examined superfluid turbulence (ST) in thermal counterflow systems, 
and have revealed a variety of physical phenomena [37|- Since the dynamics of quantized vortices 
are nonlinear and non-local, it has not been easy to quantitatively understand these observations 
on the basis of vortex dynamics. Schwarz clarified the picture of ST based on tangled vortices 
by numerical simulation of the quantized vortex filament model in the thermal counterflow [381 
139] ■ However, since the thermal counterflow has no analogy in conventional fluid dynamics, this 
study was not helpful in clarifying the relationship between ST and classical turbulence (CT). 
Superfluid turbulence is often called quantum turbulence (QT), which emphasizes the belief that 
it is comprised of quantized vortices [HI [20] . 

3.2 Bose-Einstein Crowd Superfluids 

Recall from the first section, what happens if we rotate a cylindrical vessel with a classical viscous 
fluid inside. Even if the fluid is initially at rest, it starts to rotate and eventually reaches a steady 
rotation with the same rotational speed as the vessel. In that case, one can say that the system 
contains a vortex that mimics solid-body rotation. A rotation of arbitrary angular velocity can be 
sustained by a single vortex. 

However, this does not occur in a quantum superfluid. Because of quantization of circulation, 
superfluids respond to rotation, not with a single vortex, but with a lattice of quantized vortices. 
R. Feynman noted that in uniform rotation with angular velocity the curl of the superfluid 
velocity is the circulation per unit area, and since the curl is 217, a lattice of quantized vortices 
with number density no = curl Vx/n = 2^1/ hl (the 'Feynman rule') arranges itself parallel to the 
rotation axis [29 . Such experiments were performed for superfluid ^He: Packard et al. visualized 
vortex lattices on the rotational axis by trapping electrons along the cores [40l [41]. This idea 
has also been applied to atomic Bose-Einstein condensates. Several groups have observed vortex 
lattices in rotating BECf[l!l 43 [S] [45l [42] . 

This observation has been reproduced by a simulation of the Gross-Pitaevskii (GP) equation 
for the macroscopic wave-function ^p{K,t) — \ilj{x,t)\e^^^^'*'> in 2D [46] [47] and 3D [48] spaces. 
In a weakly interacting Bose system, the macroscopic wave- function Tp(yi,t) appears as the order 
parameter of the Bose-Einstein condensate (representing our quantum crowd superfluid model) 

driven by an injected heat current, and the normal fluid and superfluid flow in opposite directions. The superflow 
was found to become dissipative when the relative velocity between the two fluids exceeds a critical value I26| . 
Gorter and Mellink attributed the dissipation to mutual friction between two fluids, and considered the possibility 
of superfluid turbulence. Feynman proposed a turbulent superfluid state consisting of a tangle of quantized vortices 
|29) . Hall and Vinen performed the experiments of second sound attenuation in rotating *He, and found that the 
mutual friction arises from interaction between the normal fluid and quantized vortices 1301 1311 ; second sound refers 
to entropy wave in which superfluid and normal fluid oscillate oppositely, and its propagation and attenuation give 
the information of the vortex density in the fluid. 

^^Among them, Madison et al. directly observed nonlinear processes such as vortex nucleation and lattice for- 
mation in a rotating ^'^Rb BEG |45l . By sudden application of a rotation along the trapping potential, an initially 
axi-symmetric condensate undergoes a collective quadrupole oscillation to an elliptically deformed condensate. This 
oscillation continues for a few hundred milliseconds with gradually decreasing amplitude. Then the axial symmetry 
of the condensate is recovered and vortices enter the condensate through its surface, eventually settling into a lattice 
configuration. 
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obeying the GP equation, or the modified cubic NLS equation (extended by the Unear term —fiip): 

indt^p^-—Aip + v\ij\'^^P-fii;. (10) 

2m 

Writing ■0 = \tp\ exp(i0), the squared amphtude is the crowd superfluid density and the gradient 
of the phase 9 gives the crowd superfluid velocity "Vx — {h/m)\J9, corresponding to frictionless flow 
of the crowd. This relation causes quantized vortices to appear with quantized crowd circulation. 
The only characteristic scale of the GP model is the coherence length defined by ^ = h/{\/2mV 
which determines the crowd vortex core size. The GP model can explain not only the crowd vortex 
dynamics but also phenomena related to vortex cores, such as crowd reconncction and nucleation. 




Figure 4: A typical crowd vortex lattice formation (modified and adapted from [T5]). 

A typical 2D numerical simulation of equation (ITOl) (adapted from [46l |47] ) for the crowd 
vortex lattice formation is shown in Figure lU where the crowd superfluid density and the phase 
are displayed together. The trapping potential is 

l/ex = \rnuj'^[{l + + (1 + ^y)v\ 

where a; = 27r x 219 Hz, and the parameters ex and describe small deviations from axisym- 
metry corresponding to experiments [HI SS] . Following [HJ [5D] , we first prepare an equilibrium 
condensate trapped in a stationary potential; the size of the condensate cloud is determined by 
the Thomas-Fermi radius i?TF- When we apply a rotation with £7 = 0.7w, the condensate becomes 
elliptic and performs a quadrupole oscillation [Fig. lUJa)]. Then, the boundary surface of the con- 
densate becomes unstable and generates ripples that propagate along the surface [Fig. IH^b)]. As 
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stated previously, it is possible to identify quantized vortices in the phase profile also. As soon as 
the rotation starts, many vortices appear in the low-density region outside of the condensate [Fig. 
m^a)]. Since quantized vortices are excitations, their nucleation increases the energy of the system. 
Because of the low density in the outskirts of the condensate, however, their nucleation contributes 
little to the energy and angular momentum!^ Since these vortices outside of the condensate are 
not observed in the density profile, they are called 'ghost vortices'. Their movement toward the 
Thomas- Fermi surface excites ripples [Fig. lUJb)]. It is not easy for these ghost vortices to enter 
the condensate, because that would increase both the energy and angular momentum. Only some 
vortices enter the condensate cloud to become 'real vortices' wearing the usual density profile of 
quantized vortices [Fig. |l{d)], eventually forming a vortex lattice [Fig. UJ^e) and (f)]. The number 
of vortices forming a lattice is given by 'Feynman's rule' no = 2U,/ n. The numerical results agree 
quantitatively with these observations. Here we remark on the essence of nonlinear dynamics. The 
initial state has no vortices in the absence of rotation. The final state is a vortex lattice corre- 
sponding to rotational frequency fi. In order to go from the initial to the final state, the system 
makes use of as many excitations as possible, such as vortices, quadruple oscillation, and surface 
waves. These experimental and theoretical results demonstrate typical behavior of quantum fluid 
dynamics in atomic BECs (THl [20] ■ 

3.3 Kolmogorov Energy Spectra 

The Kolmogorov energy spectra were confirmed for both decaying [54] and steady [55^ QT by the 
GP model. The normalized GP equation is: 

i9tV = -^Ai^ + y|V'|V-MV', (11) 

which determines the dynamics of the macroscopic wave-function ■(/'(x, t) — /(x, exp[i(/)(x, t)]. 
The crowd supcrfluid dynamics in the GP model are compressible. The total number of crowd 
agents is = / [V'l^dx and the total crowd energy is: 

as represented by the sum of the interaction energy Eint{t), the quantum energy Eq{t), and the 
kinetic energy Ekin {t) [56l [57] , 

E^nt{t) = fd^^ Eq{t) = 1 y [V/]2dx, Eu^n{t) ^^J [/V^J^dx. 

The kinetic energy is further divided into a compressible part -Efej„(t) due to compressible excita- 
tions and an incompressible part El^^{t) due to vortices. The Kolmogorov spectrum is expected 

iorEi^jtSmm- 

Actually the vortex-antivortex pairs are nucleated in the low-density region. Then the vortices parallel to the 
rotation are dragged into the Thomas-Fcrmi surface, while the antivortices are repelled to the outskirts. 

■^'^The failure to obtain a Kolmogorov law in the pure GP model 56, 57 is attributable to the following (see 
|19l I20|V The simulations showed that E^^^{t) decreases and E^^^{t) increases while the total energy E{t) is 
conserved because many compressible excitations are created through vortex reconnections 1521 1531 and disturb the 
Richardson cascade of quantized vortices. Kobayashi et al. overcame these difficulties and obtained a Kolmogorov 



16 



4 A Variety of Crowd Waves 



4.1 Crowd Shock- Waves, Solitons and Rogue Waves 

The general crowd NLS equation ([ij was exactly solved in [60l |61] [58] using the power series 
expansion method of Jacobi elliptic functions |62j . Consider the -0— function describing a single 
plane wave, with the wave number k and circular frequency w: 

i/;(a;,t) = 0(^)e'(''"^""*\ with C a; - fci and e R. (12) 

Its substitution into the NLS equation ([IJ gives the nonlinear crowd oscillator ODE: 

+ [uj- h^] m - V(^'iO = 0. (13) 

We can seek a solution 0(^) for ([T^ as a linear function [ST] 

'/'(O = flo + aisn(^), 

where sn(x) = sn(a;, m) are Jacobi elliptic sine functions with elliptic modulus m £ [0, 1], such that 
sn(a;, 0) — sin(x) and sn(x, 1) — tanh(a;). The solution of ([T3|) was calculated in 58 to be 

for m e [0, 1]; and 
for m — I. 

This gives the exact periodic solution of ([1]) as [SS] 



/ct)e'['=^-5*(i+™'+fe')l, for me [0,1); (14) 



fci)e'['="-5*(2+fe=)l, for m = l, (15) 

where (jl4[) defines the general solution, while (1151) defines the crowd envelope shock-wav^^ (or, 
'dark soliton') solution of the crowd NLS equation ([iJ. 

spectrum in QT that revealed an energy cascade 1541 1551 . By performing numerical calculations of the Fourier- 
transformed GP equation with dissipation, they confirmed the Kolmogorov spectra for decaying turbulence I54| . 
To obtain a turbulent state, they started the calculation from an initial configuration in which the density was 
uniform and the phase of the wave-function had a random spatial distribution. The initial state was dynamically 
unstable and soon developed turbulence with many vortex loops. The spectrum E''i^^^{k,t) was then found to obey 
the Kolmogorov law. 

^^A shock wave is a type of fast-propagating nonlinear disturbance that carries energy and can propagate through 
a medium (or, field). It is characterized by an abrupt, nearly discontinuous change in the characteristics of the 
medium. The energy of a shock wave dissipates relatively quickly with distance and its entropy increases. On the 
other hand, a soliton is a self-reinforcing nonlinear solitary wave packet that maintains its shape while it travels at 
constant speed. It is caused by a cancelation of nonlinear and dispersive effects in the medium (or, field). 




1pi{x,t) = ±TO. 



V{w\ 



■ sn(x 



V{w) 



tanh(a; 
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Alternatively, if we seek a solution ^(^) as a linear function of Jacobi elliptic cosine functions, 
such that cn(2:, 0) — cos(a;) and cn(x, 1) = sechfalF^ 



HO = aa + aicn(^). 



then we get 



ij^{x,t) =FmJ-j^p^cn(a;-H)e'['=^-^*(i-2™'+'=')l, for me [0,1); (16) 



'iP^{x,t) = zp /--l-sech(a;~/ct)e'['=^-^*('='-^)], formal, (17) 

where p6p defines the general solution, while (|17p defines the crowd envelope solitary-wave (or, 
'bright soliton') solution of the crowd NLS equation ([1]). 

In aU four solution expressions p^ . ([T5|). and pTl) . the adaptive crowd potential V{w) is yet 
to be calculated using either unsupervised Hebbian learning, or supervised Levenberg-Marquardt 
algorithm (see, e.g. |11)). In this way, the NLS equation ([ij becomes a quantum neural network 
[10) . Any kind of numerical analysis can be easily performed using above closed- form crowd 
solutions ip^{x,t) {i — 1, ...,4) as initial conditions. 

In addition, two new wave-solutions of the crowd NLS equation ([1]) have been recently provided 
in |63j , in the form of rogue waves using the deformed Darboux transformation method developed 
in l66l: 



1. The one-rogon crowd solution: 



^lrogon(a;,i) = « 



4(l-f-a^t) 
l + 2a2(x-fct)2 + a2a4i2 

where a and k denote the crowd scaling and gauge. 

2. The two-rogon crowd solution: 

^^P2{x,t)+lQ2{x,t) 



[kx + l/2(a'^-k'^)t] y > (18) 



e 



i[kx+l/2{a^-k^)i\^ 1/ > 0, (19) 



R2{x,t) 

where P2, Q2, R2 are certain polynomial functions of x and t. 

Both rogon crowd solutions can be easily made adaptive by introducing a set of 'synaptic 
weights' for nonlinear data fitting, in the same way as before. 



^^A closely related solution of an anharmonic oscillator ODE: 

<li"{s) + <l,{s) + <ly^{s) = 

is given by 



^, . 2m / / 2m 

= \ cn 1 / 1 H s, m 

V 1 - 2m 1 - 2m 

Rogue waves are also known as freak waves, monster waves, killer waves, giant waves and extreme waves. They 
are found in various media, including optical fibers 1641 . The basic rogue wave solution was first presented by 
Peregrine [65] to describe the phenomenon known as Peregrine soliton (or Peregrine breather). 
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4.2 Collision of Two Crowds 



Next, a bidirectional quantum neural network resembling the strong crowd coupling model ([3]) has 
been formulated in [SS] as a self-organized system of two coupled NLS equations: 

Red NLS: idta = -^9,,^ + + a, (20) 

Blue NLS: idt^p = -^9,,^ + V^(w') (Ifrp + IV'I') ^- (21) 

In this coupled model, the cr-NLS (1^(71) governs the (x, t)— evolution of the red crowd, which plays 
the role of a nonlinear coefhcient in the blue crowd (HU; the ?/>-NLS (PT|) defines the (x, i)— evolution 
of the blue crowd, which plays the role of a nonlinear coefficient in the red crowd (PO)) . The purpose 
of this coupling is to generate the crowd leverage effect (similar to stock leverage effect in which 
stock volatility is (negatively) correlated to stock returns. This bidirectional associative memory 
effectively performs quantum neural computation |10) . by giving a spatiotemporal and quantum 
generalization of Kosko's BAM family of neural networks [571 IBH]- In addition, the shock- wave 
and solitary-wave nature of the coupled NLS equations may describe brain-like effects frequently 
occurring in crowd dynamics: propagation, reflection and collision of shock and solitary waves (see 
[59]). 

The coupled crowd NLS-system ([20| -(|2T |) . without an embedded w— learning (i.e., for con- 
stant V), actually defines the well-known Manakov system^^ proven by S. Manakov in 1973 [69] 
to be completely integrable, by the existence of infinite number of involutive integrals of mo- 
tion. Manakov's own method was based on the Lax pair representation^^ It admits both 'bright' and 
'dark' soliton solutions. The simplest solution of (l20l) - ([2T]) . the so-called Manakov bright 2-soliton 
(see FigureO, has the form resembling that of the sech-solution (H?)) (see [7^ [751 [731 [7^ 1751 [771 [75] ) . 
and is formally defined by: 

V'xoi(^,0 = 2&csech[26(a; + 4at)]c-2i(2a^t+ax-26^t)^ ^23) 
where rp^^i^x^t) — I , ) \i 1 1 c = (ci,C2)^ is a unit vector such that |cip + |c2p ~ 1. Real- 



ip{x,t) 

valued parameters a and b are some simple functions of (V,fc), which can be determined by the 
Levenberg-Marquardt algorithm. 



Manakov system has been used to describe the interaction between wave packets in dispersive conservative 
media, and also the interaction between orthogonally polarized components in nonlinear optical fibres (see, e.g. 
|70l I71| and references therein) . 
^®The Manakov system lj20|l - lj2Hl has the following Lax pair |79) representation: 

dx(j> = M(j) and dt4> = B4>, or d^B - dtM = [M, B], with (22) 

/ -iA Vi ^2 \ 
M(A) = -V'l iA and 

V -^2 iA y 

2A2 - IV^iP - 2ilPi\-d^tpi 2it/>2A - 9a:l/'2 

B{X) = -i\ -2ii^lX - da^tp-* -2A2-h|,/;^|2 

-2iV'2A - a^i/'l Vi'/'2 -2A2 4-|V'2p 
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Figure 5: Hypothetical crowd-collision scenario of the Manakov 2-soliton ([23l) . Due to symmetry 
of the Manakov system, the two crowds [ip and a) can exchange their roles. 

4.3 Quantum Linear Crowd Waves 

In the case of very weak crowd heat potential V{w) <^ 1, we have V^(V') 0, and therefore equation 
([1]) can be approximated by a quantum-like crowd wave packet. It is defined by a continuous 
superposition of de Broglie's plane waves, 'physically' associated with a free quantum particle of 
unit mass. This linear wave packet, given by the time-dependent complex-valued wave function 
ip = tp{x,t), is a solution of the linear Schrddinger equation with zero potential energy and the 
crowd Hamiltonian operator H. This equation can be written as: 

idti^^H^, where H ^ ^^8^,. (24) 

Thus, we consider the i/'" function describing a single de Broglie's plane wave, with the wave 
number fc, linear momentum p = k, wavelength Afe — 27r/fc, angular frequency uik = A:^/2, and 
oscillation period Tj. — 2n/ujk = 47r/fc^. It is defined by (compare with [501 15T1 

■ipkix,t) = ^e'^'^^-"''*) = ^e'^''^"^*) =Acos(fca;- yt) + ylisin(fca;- yi), (25) 

where A is the amplitude of the wave, the angle {kx — ujkt) ~ (kx ~ ^t) represents the phase of 
the wave ip). with the crowd phase velocity: = uJk/k = k/2. 

The space-time wave function ip{x,t) that satisfies the linear Schrodinger equation can be 
decomposed (using Fourier's separation of variables) into the spatial part 4>{x) and the temporal 
part e~'"* as: 

^(x,t) = (i,{x)e-''^' = 0(x)e-'^* = 0(a;)e-5fc'*, 

where Planck's energy quantum of the crowd wave ipj. is given by: Ek = uJk = -^fc^. 

The spatial part, representing stationary (or, amplitude) wave function, (j>{x) — Ae'*'^, satisfies 
the crowd harmonic oscillator, which can be formulated in several equivalent forms: 

0" + k^(f> = 0, 0" + (^^^ = 0, (/)" + 2Ek(f> = 0. (26) 
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From the plane- wave expressions (pS)) we have: il}^{x^t) = Ae'*^^^~'^'=*^— for the wave going to 
the 'right' and ^i.{x,t) = ^e"'^*'^'*"^''*''— for the wave going to the 'left'. 

The general solution to is formulated as a linear combination of de Broglie's planar waves 
(|25p . comprising the crowd wave-packet: 

n 

ijj{x,t) ='^c^■^p^.^{x,t), (withneN). (27) 

i=0 

Its absolute square, \'ip{x,t)\'^ , represents crowd's probability density function at a time t. 

The crowd group velocity is given by: Vg = duk/dk. It is related to the crowd phase velocity 
Vk'- Vg = Vk — Xkdvk/dXk- Closely related is the center of the crowd wave-packet (the point of 
maximum crowd amplitude), given by: x — tduj^/dk. 

The following quantum-motivated assertions can be stated: 

1 . The total energy E of an crowd wave-packet is (in the case of similar plane waves) given by 
Planck's superposition of the energies Ek of n individual agents' waves: E = ntok = f fc^, 
where L — n denotes the angular momentum of the crowd wave-packet, representing the 
shift between its growth and decay, and vice versa. 

2. The average energy {E) of an crowd wave-packet is given by Boltzmann's partition function: 

Er^o nEkQ--^ _ Ek 



EDO _ 
„=o e efT -1 

where b is the Boltzmann-like kinetic constant and T is the crowd 'temperature'. 

3. The energy form of the Schrodinger equation (|24p reads: E4> = i9tV'- 

4. The eigenvalue equation for the crowd Hamiltonian operator H is the stationary Schrodinger 
equation: 

H(t){x) = E(t){x), or E4){x) ^ --d^^cpix), 

which is just another form of the harmonic oscillator (|26l) . It has oscillatory solutions of the 
form: 

Mx) = cie'^- + C2e-'^- , 
called energy eigen-states with energies Ek and denoted by: H(j3^{x) — Ek(j)E{x). 

Now, given some initial crowd wave function, -0(2;, 0) — 'iIjq{x), a solution to the initial-value 
problem for the linear Schrodinger equation (|24p is, in terms of the pair of Fourier transforms 
(J", J"-i), given by (see [SI]) 

V'(x, t) = [e-'-* J-(V^o)] = [e-^'^'Ti^Po)] ■ (28) 

For example (see |81j). suppose we have an initial crowd wave- function at time t — given by 
the complex-valued Gaussian function: 

V'(x,0) =e-'"^'/2gifc:.^ 



21 



where a is the width of the Gaussian, while p is the average momentum of the wave. Its Fourier 
transform, ■0o(fc) = J"[^/;(a;, 0)], is given by 



(fc-p)^ 



The sohition at time t of the initial value problem is given by 



'2'Ka 

which, after some algebra becomes 

Vl + iai 

As a simpler example!^ if we have an initial crowd wave-function given by the real-valued 
Gaussian function. 



^(a:,0) 



the solution of (l24l) is given by the complex-valued V'— function, 

,/ +^ °^P(~2(i+it)) 

From (j28p it follows that a stationary crowd wave-packet is given by: 
1 f^°° 

(t){x) = —= \ e'*^"' where = 

V 27r J-oc 



As is the stationary crowd PDF, we can calculate the crowd expectation values and the 

wave number of the whole crowd wave-packet, consisting of n measured plane waves, as: 

x\^{x)\'^dx and (fc) = / k\ij{k)\'^dk. (29) 

00 J —oc 

The recordings of n individual crowd plane waves (j25p will be scattered around the mean values 
(|29p. The width of the distribution of the recorded x— and fc— values are uncertainties Aa; and 
AA;, respectively. They satisfy the Heisenberg-type uncertainty relation: 

AxAk > -, 
- 2' 
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An example of a more general Gaussian wave-packet solution of H24| l is given by: 



ip(x,t) = \ — — : — exp ' 



1 + iat \ 1 -f iat 

where soiPO a-re initial stock- price and average momentum, while a is the width of the Gaussian. At time t = 
the 'particle' is at rest around s = 0, its average momentum po = 0. The wave function spreads with time while 
its maximum decreases and stays put at the origin. At time —t the wave packet is the complex-conjugate of the 
wave-packet at time t. 
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which imply the similar relation for the total crowd energy and time: 

71 

AEAt> -. 
- 2 

5 Conclusion 

In this paper we gave a formal mathematical and physical description of nonlinear phenomena 

in dynamics of human crowds. While Helbing discovered a phenomenon of crowd turbulence 
(see Introduction), we felt that equally important would be to model related but different crowd 
phenomena, such as solitons, rogue waves and shock waves. Our proposal, including both classical 
and quantum description of crowd turbulence, as well as both nonlinear and quantum crowd waves, 
provides a new basis for studying all these nonlinear phenomena in crowds. 

6 Appendix: Basic Lie Algebra Mechanics 

A manifold M is a topological space that on a small scale (locally) resembles the Euclidean space. 
Manifolds are usually endowed with a differentiable structure that allows one to do calculus and 
differential equations, as well as a Riemannian m,etric that allows one to measure distances and 
angles. For example, Riemannian manifolds are the configuration spaces for Lagrangian mechanics, 
while symplectic manifolds are the phase spaces in the Hamiltonian mechanics. A diffeomorphism 
is an invertible function that maps one smooth (differentiable) manifold to another, such that both 
the function and its inverse are smooth. 

A Lie group G is a smooth manifold M that has at the same time a group G— structure 
consistent with its manifold M— structure in the sense that group multiplication fi : G x G ^ 
G, {g, h) I— >■ gh and the group inversion v : G ^ G, g '-^ g^^ are smooth maps. A point e G G 
is called the group identity element. 

A Lie group can act on a smooth manifold M by moving the points of M, denoted by G x M — 
M. Group action on a manifold defines the orbit of a point to on a manifold M, which is the set of 
points on M to which to can be moved by the elements of a Lie group G. The orbit of a point m 
is denoted by Gto = {g ■ m\g S G}. 

Let G be a real Lie group. Its Lie algebra g is the tangent space TGe to the group G at 
the identity e provided with the Lie bracket {commutator) operation [X, y], which is bilinear, 
skew-symmetric, and satisfies the Jacobi identity (for any three vector fields X,Y, Z G q): 

[[X,Y],Z] = [X, [Y,Z]]-[X, [Y,Z]]. 

Note that in Hamiltonian mechanics, Jacobi identity is satisfied by Poisson brackets, while in 
quantum mechanics it is satisfied by operator commutators. 

For example, G = SO{3) is the group of rotations of 3D Euclidean space, i.e. the configuration 
space of a rigid body fixed at a point. A motion of the body is then described by a curve g = g{t) 
in the group SO{3). Its Lie algebra q — so{3) is the 3D vector space of angular velocities of all 
possible rotations. The commutator in this algebra is the usual vector (cross) product. 

A Lie group G acts on itself by left and right translations: every element g G G defines 
difii'eomorphisms of the group onto itself (for every h £ G): 

Lg : G ^ G, Lgh = gh; Rg : G ^ G, Rgh = hg. 
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The induced maps of the tangent spaces are denoted by: 

Z/g* : TGh —> TGgh, Rg* '■ TGh TGhg- 

The diffeomorphism Rg-iLg is an inner automorphism of the group G. It leaves the group 
identity e fixed. Its derivative at the identity e is a linear map from the Lie algebra g to itself: 

Adg : g, Adg{Rg-lLg)*e 

is called the adjoint representation of the Lie group G. 

Referring to the previous example, a rotation velocity () of the rigid body (fixed at a point) is 
a tangent vector to the Lie group G = SO{3) at the point g G G. To get the angular velocity, we 
must carry this vector to the tangent space TGg of the group at the identity, i.e. to its Lie algebra 
g = so(3). This can be done in two ways: by left and right translation, Lg and Rg. As a result, 
we obtain two different vector fields in the Lie algebra so (3) : 

LJc = Lg-i^g G so{3) and ujx = Rg-i*g e so{3), 

which are called the 'angular velocity in the body' and the 'angular velocity in space,' respectively. 

Now, left and right translations induce operators on the cotangent space T*Gg dual to Lg* and 
i?g*, denoted by (for every h £ G): 

g ■ J '^gh — >■ J ^h, ^g ■ J (j/ig ^ ^h- 

The transpose operators Ad*g : — >■ satisfy the relations Ad*gf^ = Ad\Ad*g (for every g,h G G) and 
constitute the co-adjoint representation of the Lie group G. The co-adjoint representation plays 
an important role in all questions related to (left) invariant metrics on the Lie group. According 
to A. Kirillov, the orbit of any vector field X in a Lie algebra in a co-adjoint representation Ad* 
is itself a symplectic manifold and therefore a phase space for a Hamiltonian mechanical system. 

A Ricmannian metric on a Lie group G is called left-invariant if it is preserved by all left 
translations Lg, i.e., if the derivative of left translation carries every vector to a vector of the same 
length. Similarly, a vector field X on G is called left-invariant if (for every g & G) L*gX = X. 

Again referring to the previous example of the rigid body, the dual space g* to the Lie algebra 
Q = SO (3) is the space of angular momenta tt. The kinetic energy T of a body is determined by 
the vector field of angular velocity in the body and does not depend on the position of the body 
in space. Therefore, kinetic energy gives a left-invariant Riemannian metric on the rotation group 
G = SO{3). 
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